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A  mathematical  model  was  developed  to  simulate  the  performance  of  a  direct  carbon  fuel  cell.  The  model 
takes  account  of  the  electrochemical  reaction  dynamics,  mass-transfer  and  the  electrode  processes.  An 
improved  packed  bed  anode  was  adopted.  Polarization  losses  for  the  cell  components  were  examined 
supposing  graphite  as  the  fuel  and  molten  carbonate  as  the  electrolyte.  The  results  indicated  that  the 
anode  activation  polarization  was  the  major  potential  loss  in  923-1023  K.  The  effects  of  temperature, 
anode  dimension,  and  carbon  particle  size  on  the  cell  performance  were  investigated.  The  model  predicted 
that  the  power  density  can  be  as  high  as  200-500  Wrrr2,  with  carbon  particle  size  in  the  range  1.0  x  10  7 
to  1.0  x  10-4  m  and  in  923-1023  K  and  that  the  overall  efficiency  of  the  cell  is  higher  than  55%  for  low 
current  density  and  is  45-50%  for  high  current  density. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Factors  that  impact  the  development  of  an  energy  system 
include  the  availability  of  fuel,  efficiency  and  environmental  effect 
[1  ].  Fuel  cells  have  a  high  efficiency  for  converting  chemical  energy 
to  electrical  energy.  These  devices  are  not  limited  by  theoretical 
Carnot  cycle  efficiencies  [2].  They  also  typically  do  not  emit  pollu¬ 
tants,  such  as  NOx,  SOx,  etc.  [3]. 

Since  fuel  cells  were  discovered  by  William  Grove  in  1839,  a 
family  of  fuel  cells  has  been  developed  [4],  which  include  gaseous 
hydrogen,  methane,  methanol,  ethanol,  and  solid  carbon  [5-10]. 
Flydrogen  fuel  cells  are  technically  mature  and  have  been  applied  in 
space  explorations  and  are  starting  to  replace  the  internal  combus¬ 
tion  engines  in  automobiles  [11  ].  There  has  been  renewed  interest 
on  direct  carbon  fuel  cell  (DCFC),  which  uses  solid  carbon  as  the 
fuel,  due  to  their  potential  to  use  coal  and  coal  derivatives  as  fuel 
and  high  thermal  efficiency  [8-10].  The  principle  reaction  is 

C  +  02  C02  AH973  =  395  kj  mol”1  (1 ) 

The  theoretical  thermal  efficiency  of  DCFC  is  almost  100%  [8,10]. 
The  DCFC  is  a  high  temperature  device  and  shares  features  with  a 
molten  carbonate  fuel  cell  (MCFC)  and  a  solid  oxide  fuel  cell  (SOFC) 
[12-21]. 
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The  first  DCFC  prototype  was  established  by  Jacques  more 
than  110  years  ago  [13].  Flowever,  the  research  on  this  device 
was  intermittent  in  the  early  20th  century  because  of  substan¬ 
tial  improvements  in  the  efficiencies  of  steam  power  plants.  In  the 
1970s,  the  concept  was  verified  by  Weaver  and  coworkers  [14,15]. 
Since  then,  there  has  been  substantial  progress  in  the  develop¬ 
ment  of  high  temperature  fuel  cells  including  MCFCs  and  SOFCs. 
Concepts,  materials  and  processes  proved  with  these  fuel  cells 
were  applied  to  DCFC  and  lead  to  some  major  breakthroughs  [8,9]. 
Cherepy  et  al.  [8]  investigated  the  effect  of  carbon  nanostructure  on 
the  electrochemical  reactivity  and  found  that  “turbostratic”  carbon, 
with  a  high  degree  of  disorder  in  structure  and  domain  sizes  in  the 
range  of 30-100  nm,  performed  the  best.  They  tested  a  five-cell  self¬ 
feeding  bipolar  stack  and  designed  a  larger  scale  system  for  clean 
coal.  Tao  [19]  employed  a  sophisticated  configuration  which  used  a 
liquid  tin  anode  and  a  YSZ  electrolyte.  Lipilin  et  al.  [21  ]  established 
a  device  combining  the  features  of  SOFC  and  MCFC  technology,  in 
which  a  circulating  liquid-molten-salt  electrolyte  and  a  carbon-rich 
fuel  were  used,  and  a  power  density  lOOmWcm-2  was  achieved. 
Steinberg  [20],  Flemmes  and  Flouwing  [22]  carried  out  macroscopic 
energy  and  exergy  analyses  of  a  DCFC  system.  Agarwal  and  Korn- 
hauser  [23]  examined  an  integrated-DCFC-FIPBR  ( hydrogen  plasma 
black  reactor)  system.  However,  these  modeling  studies  did  not 
account  for  the  electrode  processes  in  the  DCFC.  Understanding 
of  the  electrode  reactions,  mass  and  heat  transfers  in  the  device, 
and  the  efficiency  analyses  models  for  DCFC  including  the  electrode 
processes  are  needed. 
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Nomenclature 

A  a  kinetic  parameter  related  with  Ohmic  resistance 

(Q  m2  Pa0-67) 

b  diameter  of  the  contact  area  between  two  carbon 

particles  (m) 

B  a  kinetic  parameter  related  with  Ohmic  resistance 

(K) 

C  concentration  (mol  m-3 ) 

CCo2  C02  concentration  in  gas  phase  (mol  m-3 ) 

Cr  a  constant  related  with  Ohmic  resistance  of  the  con¬ 
tacts  (£2m2) 

D  a  kinetic  parameter  related  with  electrolyte  content 

in  the  matrix  (Q  m2) 

Dc  diameter  of  carbon  particle  (m) 

Ea  apparent  activation  energy  of  the  electrolyte  con¬ 

ductivity  (J  mol-1 ) 

/  a  kinetic  parameter  related  with  Ohmic  resistance 

(I<) 

F  mechanical  force  (N) 

F  Faraday’s  constant  (96,487  C  mol-1 ) 

Ff  buoyancy  force  (N) 

Fg  gravitational  force  (N) 

Fr  repulsive  force  (N) 

G  Gibbs-free  energy  (J) 

H  enthalpy  (J) 

i  order  number  of  a  slab 

/  current  (A) 

If  measured  current  (A) 

It  theoretical  current  (A) 

j  current  density  (A  m-2 ) 

jc  current  density  in  the  electrode  phase  (A  m-2 ) 

jlim  limiting  current  density  (A  m-2 ) 

js  current  density  in  the  electrolyte  phase  (A  m-2 ) 

jo  exchange  current  density  (A  m-2 ) 

j°  concentration-independent  exchange  current  den¬ 

sity  (Am-2) 

I<  mass  transport  coefficient  (m  s-1 ) 

n  numbers  of  electrical  charge  transfer  in  an  electro¬ 

chemical  reaction  (mole) 

N\  number  of  the  carbon  particles  in  each  slab 

Ns  number  of  slabs 

pm  mean  yield  pressure  of  the  asperities  (Pa) 

P  pressure  in  gas  phase  (Pa) 

Pc o2  partial  pressure  of  C02  in  gas  phase  (Pa) 

Po2  partial  pressure  of  02  in  gas  phase  (Pa) 

r\  power  of  C02  partial  pressure  for  the  exchange  cur¬ 

rent  density  equation  in  cathode 
r2  power  of  02  partial  pressure  for  the  exchange  cur¬ 

rent  density  equation  in  cathode 
R  resistance  (£2  m2) 

Rc  resistance  in  the  electrode  phase  ( £2  m2 ) 

Fee  resistance  of  all  components  except  for  anode 
(Q  m2) 

Fcr  constriction  resistance  of  the  electrode  phase 
(£2  m2) 

Fct  charge  transfer  resistance  ( £2  m2 ) 

Ff  film  resistance  ( £2  m2 ) 

Fs  resistance  in  the  electrolyte  ( £2  m2 ) 

R'  gas  constant  (8.314  J  mol  K-1 ) 

s  contact  area  between  two  carbon  particles  (m2) 

S  entropy  (JK-1) 

T  temperature  (I<) 


Vcell  cell  voltage  (V) 

VN  open  circuit  voltage  (V) 

V°  open  circuit  voltage  at  standard  state  in  Eq.  (10)  (V) 

x  height  of  the  packed  bed  anode  (m) 

y  length  of  the  packed  bed  anode  (m) 

z  width  of  the  packed  bed  anode  (m) 

Greek  letters 

P  apparent  charge  transfer-coefficient 

y  efficiency 

Yf  current  efficiency 

Yt  theoretical  reversible  efficiency 

yv  voltage  efficiency 

e  voidage  of  the  packed  bed 

Tj  polarization  loss  (V) 

r] art  activation  polarization  (V) 

rjcon  concentration  polarization  (V) 

?7ohm  Ohmic  polarization  (V) 

cr c  conductivity  of  the  electrode  (S  m-2 ) 

crs  conductivity  of  the  electrolyte  (S  m-2 ) 

of  a  pre-exponential  factor  (S  m-2 ) 

Superscript 
b  bulk 

0  standard  state 

Subscript 

act  activation 

an  anode 

cat  cathode 

con  concentration 

lim  limit 

ohm  Ohmic 


Work  described  in  this  paper  focuses  on  the  development  of 
a  preliminary  mathematical  model  describing  the  performance  of 
a  DCFC  single  cell.  The  purpose  is  to  determine  and  simulate  the 
various  polarizations,  V-j  relation,  power  density  and  efficiencies 
of  the  cell.  A  commercial  software  MATLAB™  is  adopted  for  this 
purpose. 

2.  Model 

2.1.  System  configuration 

Previous  literature  [8]  assumed  that  the  collision  of  carbon  par¬ 
ticles  in  the  anode  transfers  the  current.  Flowever,  this  mechanism 
requires  an  unrealistically  high  collision  frequency  to  account  for 
the  experimental  results  [24].  In  this  work,  we  propose  an  alterna¬ 
tive  packed  bed  model  which  offers  a  simplified  configuration  for 
the  DCFC  anode. 

Fig.  1  gives  an  illustration  of  a  DCFC  single  cell.  It  consists  of 
three  key  components,  the  anode,  the  cathode,  and  the  electrolyte 
matrix.  Solid  carbon  fuel  is  fed  into  the  anode  in  the  form  of  particles 
and  wetted  with  molten  salt.  Each  particle  is  assumed  to  act  as  a 
rigid  sphere  and  packed  with  a  simple  hexagonal  pattern.  The  bed 
consists  of  several  slabs  of  spherical  carbon  particles.  The  number 
Ns  of  the  slabs  is  dependent  on  the  bed  height  x  and  the  diameter 
of  the  carbon  particle  D( 


(2) 
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Fig.  1.  (a)  Schematic  configuration  of  the  DCFC  single  cell  using  packed  bed  anode,  (b)  top  view  of  packed  bed  electrode. 


The  number  N\  of  the  particles  in  each  slab  is 

N,=  ld^ 

D2C 


(3) 


The  anode  reaction  takes  place  at  the  interface  of  carbon  and 
carbonate  melt,  with  release  of  C02  and  electrons.  The  electrolyte 
matrix  is  made  of  a  porous  ceramic  membrane  filled  with  molten 
carbonate,  which  allows  the  migration  of  carbonate  ions  from  the 
cathode  compartment  to  the  anode  compartment.  The  structure  of 
the  cathode  is  similar  to  the  one  adopted  in  a  MCFC.  A  mixture  of  air 
and  C02  with  a  ratio  of  02  to  C02  1 :2  is  introduced  into  the  cathode. 


2.2.  Electrode  processes 


The  C032-  ions  travel  across  the  electrolyte  matrix  to  the  anode  and 
the  overall  cell  reactions  are 

C  +  02  -*  C02  (8) 

2C  +  02  -*  2 CO  (9) 

The  ratio  of  C02/C0  produced  depends  on  the  thermodynamic 
and  kinetic  factors.  According  to  published  experimental  results, 
C02  is  the  major  gas  released  from  anode  in  the  temperature  range 
923-1073  K  and  the  C02  concentration  exceeds  that  predicted  by 
thermodynamic  equilibrium  [25-27].  In  this  work,  the  production 
of  CO  in  the  anode  is  not  considered.  Therefore,  reactions  (5),  (6) 
and  (9)  are  ignored  and  the  anode  reactions  are  simplified  as  a  four- 
electron  (n  =  4)  process,  i.e.  reaction  (1). 


Molten  carbonate  is  a  promising  electrolyte  for  the  electrochem¬ 
ical  oxidation  of  carbon  [8,14,15].  In  this  work,  a  Li2C03/K2C03 
eutectic  melt  with  a  molar  ratio  of  32:68  is  employed  as  the  elec¬ 
trolyte.  Graphite  is  employed  as  the  anode  fuel.  The  anode  reactions 


are  as  follows: 

C  +  2C032-  3C02  +4e-  (4) 

2C  +  C032-  3CO  +  2e~  (5) 

The  reverse  of  Boudouard  reaction 

C02  +  C^2C0  (6) 

is  a  major  side  reaction  and  lowers  the  efficiency.  C032-  ions  are 
produced  in  the  lithium  doped  NiO-based  porous  cathode  where 
the  following  reduction  reaction  occurs 

2C02+02+4e-^  2C032-  (7) 


2.3.  Electrochemical  thermodynamic  equations 

The  open  circuit  voltage  VN  is  calculated  with  Nernst  equation 
that  is  given  by 


1/  i/O  R'T , 

V  +  ^7  ln 


(  PQ2PcQ2,cat\ 
\  Pc02,an  ) 


(10) 


The  open  circuit  voltage  at  standard  state  V°  can  be  calculated  by 


Vu  = 


AG° 

h W 


(11) 


2.4.  Polarization  losses 


When  current  flows  in  the  cell,  the  voltage  decreases.  This  phe¬ 
nomenon  is  well  known  as  the  polarization  loss  q,  which  is  caused 
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by  different  mechanisms  including  charge  transfer  pact,  ion  concen¬ 
tration  gradient  ijcon ,  and  circuit  Ohmic  resistance  r]ohm.  Therefore, 
the  cell  voltage  r]ce\\  is  determined  by  subtracting  these  polarization 
losses  from  the  open  circuit  voltage 

^cell  —  —  (?7act,an  +  ?7act,cat  +  Pohm  +  ^7con,an  +  ??con,cat)  (12) 


2.4.1  Activation  polarization 

In  order  to  simulate  the  charge  transfer  processes,  the  general 
Bulter-Volmer  equation  was  used  to  calculate  the  activation  polar¬ 
izations  for  anode  and  cathode,  respectively 


j  =jo  |  exp  (ftnJ^act)  -  exp 


-(1-/8) 


nJFr}} 


(13) 


The  apparent  charge  transfer-coefficient  /3  takes  value  in  the  range 
0-1,  depending  on  the  symmetry  of  the  transition  state  in  the  elec¬ 
trochemical  reaction,  and  is  typically  assigned  as  0.5  for  fuel  cells 
[28].  The  anode  exchange  current  density  j0,an  is  related  to  the 
charge  transfer  resistance  Rc t  and  is  given  by 


jo,  an 


R'T 

nF'Rct 


(14) 


and 


Ret  =  30.79  +  1.48  exp  -gg)  (15) 

The  values  in  Eq.  (15)  are  obtained  from  a  regression  analysis  on 
the  data  published  by  Peelen  et  al.  [29]. 

For  the  cathode,  there  is  no  model  that  fits  all  the  experimental 
data.  In  this  work,  the  model  based  on  the  peroxide  mechanism  is 
selected  [30] 


Jo,cat=jo,cat(Pc02)r’(Po2)r2 


(16) 


The  values  of  cat,  V\  and  r2  were  taken  from  [31  ]. 


The  contact  area  s  is  dependent  on  the  force  F  acting  on  the  particles 
and  the  mean  yield  pressure  of  the  asperities  pm 


st  =  nb\ 


(22) 

(23) 


The  forces  applied  on  the  particles  include  gravitational  force 
Fg)j,  buoyancy  force  Ff  ?-  and  repulsive  force  Fr  i.  However,  the  repul¬ 
sive  force  on  the  top  and  the  bottom  of  carbon  particle  is  assumed 
to  be  equal.  Consequently,  the  difference  of  gravitational  force  and 
buoyancy  force  is  the  force  pressing  the  particles 


Fj  ~  (Ns  -  1  )(Ffj;  -  Fgj) 

(24) 

FfJ  =  1.33psjrg(Dc)2 

(25) 

Fgj  =  \.33pc7rg(Dc)2 

(26) 

The  resistances  of  the  other  parts  can  be  expressed  as  follows 
[35]: 

R«  =  Ayexp^+CR  +  D  exp(f)  (27) 

Fo2 

where  po2  is  the  partial  pressure  of  oxygen  at  cathode  surface,  Cr 
is  a  constant  related  with  Ohmic  resistance  of  the  contacts,  D  is  a 
kinetic  parameter  which  is  proportional  to  the  electrolyte  content 
of  the  matrix.  The  overall  Ohmic  polarization  can  be  express  as 


*7ohm  —  jR  ce  +  Pohm,an  (28) 

Ns  Ns 

*7ohm,an  —  ^  Js,jRs,i  +  ^  Jc,jRc,i  (29) 

i=l  i=l 


2.4.2.  Ohmic  polarization 

Ohmic  polarization  is  caused  by  resistances  inside  the  cell  and 
obeys  Ohm’s  law.  For  the  anode,  the  resistance  in  each  slab  con¬ 
tains  two  parts,  i.e.  the  resistances  in  the  electrolyte  Rs  i  and  in  the 
electrode  RCti  in  slab  i,  respectively  [32].  The  Rs<i  is  expressed  as 


Rs,i  = 


orsyz£l 


(17) 


The  conductivity  of  the  electrolyte  crs  is  a  function  of  tem¬ 
perature  and  obeys  Arrhenius’  law.  Here  the  parameters  were 
extrapolated  to  923  K  [33] 


ffs  =  (J“exp  (-Jy) 


(18) 


2.4.3.  Concentration  polarization 

The  reactant  transport  resistances  lead  to  a  concentration  polar¬ 
ization.  The  concentration  polarization  can  be  measured  using  a 
limiting  current  density  method  as  described  below  [36] 


In  a  DCFC,  the  phenomenon  of  current  saturation  does  not  occur, 
because  solid  carbon  is  always  in  excess  at  the  anode  [26,29].  As  a 
consequence,  the  concentration  polarization  of  the  anode  is  negli¬ 
gible.  The  cathode  is  porous,  and  the  partial  pressures  of  gases  at 
the  interface  are  lower  than  in  the  bulk  gas  phase  when  the  cur¬ 
rent  density  is  high.  The  cathode  concentration  polarization  can  be 
calculated  as 


The  current  flow  through  the  electrolyte  phase  is 


Js,i  —  J 


1-2  71 


fe)] 


Jlim,cat  —  ftF  ^C02^cq2  (31) 

The  mass  transport  coefficient  I<Co2  is  taken  as  3.5  x  10-2  ms-1  [37]. 


and  the  current  flow  through  the  electrode  phase  is 

=  2jTj  (d~t)  (20) 

Rc  i  consists  of  the  constriction  resistance,  Rcr<i,  and  the  film  resis¬ 
tance,  Rf  i  on  the  electrode  surface.  However,  for  the  fresh  carbon 
particle,  the  contact  resistance  of  two  carbon  particles  is  the  same  as 
its  constriction  resistance  and  is  related  to  the  conductivity  of  the 
carbon  particle  crc  and  the  diameter  of  the  contact  area  between 
two  particles  b  [34] 


Rc,i  —  R 


cr  ,i 


1 

orc  hi 


(21) 


2.5.  Efficiencies 


The  operating  efficiency  is  lowered  by  activation,  Ohmic,  and 
mass  transport  polarizations.  The  theoretical  reversible  efficiency 
of  a  fuel  cell  can  be  expressed  as  the  function  of  the  total  free  energy 
and  enthalpy  changes: 


Yt- 


AG 

AH 


The  overall  efficiency  is  defined  as: 


(32) 


Y  =  YtYvYf 


(33) 
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Table  1 

Setting  parameters  used  in  simulation 


Parameter 

Value 

Reference 

Total  gas  pressure  in  anode  (Pa) 

1.01325  x  105 

Total  gas  pressure  in  cathode  (Pa) 

1.01325  x  105 

Height  of  packed  bed  anode,  x  (m) 

0 

X 

q 

0 

r" 

O 

X 

q 

Length  of  packed  bed  anode,  y  (m) 

1.0  x  10-2 

Width  of  packed  bed  anode,  z  (m) 

1.0  x  10~2 

Diameter  of  spherical  graphite  particle, 

1.0  x  10-7  to  1.0  x  10-4 

Dc  (m) 

Voidage  of  packed  bed  anode,  e 

0.396 

Cathode  material 

LixNi!_xO 

[30] 

Cathode  thickness  (m) 

0.5  x  10-3 

[30] 

Electrolyte  composition  (mol%) 

62U2CO3/38K2CO3 

Matrix  thickness  (m) 

0.5  x  10"3 

[30] 

Matrix  material 

7-UAIO2 

[30] 

Cell  operating  temperature  T  (I<) 

923-1023 

Fuel 

Graphite 

Conductivity  of  graphite  particle,  crc 

1.6  x  105 

[36] 

(sm-2) 

Density  of  graphite  particle,  pc  (kgm-3) 

1.5  x  103 

[36] 

Density  of  electrolyte,  ps  (kgm-3) 

2.0  x  103 

[33] 

Mass  transport  coefficient  of  CO2,  I<co2 

3.5  x  10~2 

[32] 

(ms-1) 

Apparent  activation  energy  Ea  for  Eq. 

2.191  x  104 

[33] 

(18)  (J  mol-1) 

Pre-exponential  factor  of  for  Eq.  (18) 

2.322  x  105 

[33] 

(sirr2) 

Apparent  charge  transfer-coefficient,  ft 

0.5 

Concentration-independent  exchange 

5.0  x  102 

[35] 

current  density,  (Am-2) 

Power  of  CO2  partial  pressure  for  the 

-1.250 

[34] 

exchange  current  density  equation  in 
cathode,  r\ 

Power  of  O2  partial  pressure  for  the 

0.375 

[34] 

exchange  current  density  equation  in 
cathode,  r2 

where  yv  is  the  voltage  efficiency  given  by: 


and  yf  is  the  current  efficiency  written  as: 

*-'i  <35> 

2.6.  Further  assumptions 

All  parameters  used  in  this  work  are  summarized  in  Table  1. 
Three  additional  assumptions  for  the  simulation  include: 

•  Air  consists  of  75.67%  N2,  20.35%  02,  3.03%  H20,  0.92%  Ar,  and 
0.03%  C02  and  all  the  components  are  regarded  as  ideal  gas. 

•  The  operation  temperature  of  the  DCFC  is  in  the  range 
923-1023  K. 

•  The  whole  surface  of  the  spherical  graphite  particle  is  effective  in 
the  anode. 

3.  Results  and  discussion 

3.1.  Thermodynamic  properties 

The  open  circuit  voltage  VN  is  a  primary  parameter  for  the  per¬ 
formance  of  a  fuel  cell.  Generally,  DCFC  operates  at  the  atmospheric 
pressure  and  air  is  used  in  the  cathode.  Consequently,  the  open  cir¬ 
cuit  voltage  Vn  depends  only  on  the  working  temperature.  As  shown 
in  Fig.  2,  VN  decreases  slightly  with  the  increase  of  the  temperature 
and  is  beyond  0.87  V  in  873-1173  K. 


Fig.  2.  Open  circuit  voltage  vs.  temperature. 


3.2.  Electrochemical  dynamic  characteristics 

Figs.  3-6  illustrate  the  contributions  of  various  polarizations  as  a 
function  of  temperature  for  a  single  cell  in  923-1023  K  and  based  on 
x  =  1.0  x  10-3  m,y  =  z=  1.0  x  10-2  m,  Dc  =  1.0  x  10-5  m.  These  results 
indicate  that  the  cathode  activation  polarization  loss  is  close  to 
the  cathode  concentration  polarization,  and  both  are  much  lower 
than  the  anode  activation  polarization  loss.  The  anode  Ohmic  polar¬ 
ization  is  in  the  same  order  of  magnitude  as  the  anode  activation 
polarization. 

3.3.  Effect  of  temperature 

The  open  circuit  voltage  decreases  about  0.015  V  along  with 
the  temperature  increases  from  923  to  1023  K.  All  the  polarization 
losses  examined  show  strong  dependence  on  the  temperature.  The 
increase  of  temperature  lowers  the  anode  activation  loss  effectively, 
but  increases  the  cathode  activation  loss  slightly.  The  concentra¬ 
tion  loss  increases  a  little  with  the  increase  of  the  temperature, 
while  the  Ohmic  polarization  loss  shows  a  slight  dependence  on 
the  temperature. 


Fig.  3.  Anode  activation  polarization  loss  vs.  current  density  in  923-1023  K.  )  at 
923  K,  (O)  at  973  K,  (>)  at  1023  K. 
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Fig.  4.  Cathode  activation  polarization  loss  vs.  current  density  in  923-1023  K.  (★  ) 
at  923  K,  (O)  at  973  K,  (>)  at  1023  K. 


Fig.  5.  Cathode  concentration  polarization  loss  vs.  current  density  in  923-1023  K. 
(* )  at  923  K,  (O)  at  973  K,  (>)  at  1023  K. 


Fig.  6.  Ohmic  polarization  loss  vs.  current  density  in  923-1023 1<.  (★ )  at  923  K,  (O) 
at  973  K,  (>)  at  1023  K. 
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Fig.  7.  Effect  of  height  of  packed  bed  anode  on  Ohmic  polarization  in  slab  Nt  at  973  K 
with  Dc  =  1  x  10“5  m.  (a)  Ns  =  10;  (b)  Ns  =  50;  (c)  Ns  =  100. 


3.4.  Effect  of  anode  dimension  and  particle  size 

The  voltage  drop  in  the  anode  electrode  contributes  the  major 
part  of  the  Ohmic  polarization  losses  from  all  the  components, 
which  depends  on  the  anode  dimension  and  the  carbon  particle 
size.  Here  we  examine  three  parameters  of  the  anode  dimensions: 
the  height  x ,  the  lengthy  and  the  width  z  of  the  packed  bed  anode. 
Because  each  carbon  particle  acts  like  an  electrode,  the  voltage  drop 
in  the  anode  is  independent  of  y  and  z.  As  a  consequence,  the  x 
decides  the  anode  voltage  drop.  In  Fig.  7,  the  Ohmic  polarization 
loss  in  each  slab  is  plotted  as  the  function  of  current  density  and 
bed  height  x.  The  sum  of  voltage  drop  in  each  slab  is  the  total  anode 
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Fig.  8.  Effect  of  carbon  particle  size  on  anode  Ohmic  polarization  at  973 1<.  ) 

Dc  =  1  x  10-6  m,  (O)  Dc  =  1  x  10“5  m,  (>)  Dc  =  1  x  10“4  m. 

Ohmic  polarization  in  the  anode.  As  the  bed  height  increasing,  the 
voltage  drop  in  each  slab  increases. 

The  effect  of  carbon  particle  size  on  the  anode  Ohmic  polariza¬ 
tion  is  shown  in  Fig.  8.  In  this  work,  the  range  of  carbon  particle 
size  is  1.0  x  10-7  to  1.0  x  10-4  m.  The  anode  Ohmic  polarization  is 
<0.05  V  when  the  diameter  of  carbon  particle  is  <1.0  x  10-5  m  and 
has  the  same  order  of  magnitude  as  the  Ohmic  polarization  caused 
by  other  components.  In  the  range  of  1.0  x  10-5  to  1.0  x  10-4  m,  the 
anode  ohmic  polarization  displays  a  sharp  increase  as  the  increase 
of  the  particle  size  and  contributes  the  major  part  of  Ohmic  polar¬ 
ization.  If  the  diameter  of  carbon  particle  is  >1.0  x  10-4  m,  the  anode 
polarization  will  become  much  more  serious  and  it  is  not  in  the 
scope  of  this  model. 

3.5.  V-j  relations  and  power  density 

Fig.  9  plots  the  V-j  relations  and  powder  densities  of  the 
single  cell  at  different  temperatures.  For  the  two  cases  with 
Dc<1.0xl0~4m  and  the  current  density  <200 Am-2,  the  volt¬ 
age  is  higher  than  0.5  V  in  923-1023  K.  If  the  current  density  is 
>200  Am-2,  the  voltage  is  0.4-0.5  V  for  the  cell  using  carbon  par¬ 
ticles  with  Dc  <  1.0  x  10-5  m.  Flowever,  for  Dc  >  1.0  x  10-5  m,  the 
voltage  is  lower  than  0.4  V.  For  Dc  <  1.0  x  10-5  m,  the  power  density 
in  923-1023 1<  increases  monotonically  with  the  increase  of  the  cur¬ 
rent  density,  while  for  Dc  =  1.0  x  10-4  m,  the  power  density  shows 
peaks  at  150.2,  178.0  and  202.0 Wm~2  for  923,  973  and  1023 K, 
respectively.  Also  is  shown  in  Fig.  9,  the  cell  voltage  and  the  power 
density  show  some  dependency  on  the  working  temperature. 

The  results  presented  here  indicate  that  the  performance  of  a 
DCFC  single  cell  is  under  the  control  of  the  activation  polariza¬ 
tion  and  the  Ohmic  polarization.  The  anode  activation  polarization 
accounts  for  the  most  serious  loss,  which  can  be  attributed  to  the 
low  reaction  rate  of  the  electrochemical  oxidation  of  graphite.  The 
exchange  current  density  of  anode  in  DCFC  is  below  1 0  A  m-2 ,  which 
is  rather  low  if  compared  to  that  of  the  anode  in  a  MCFC  in  the  order 
of  several  hundred  Am-2.  Utilization  of  carbon  with  good  activity 
for  electrochemical  oxidation  reaction  is  a  feasible  way  to  overcome 
this  problem  [8]. 

3.6.  Efficiencies 

The  theoretical  reversible  efficiency  of  reaction  (1)  is  100% 
based  on  lower  heating  value  for  923,  973  and  1023  K.  While  for 
high  temperature  hydrogen  fuel  cell,  it  is  only  83%  because  of  the 


Fig.  9.  Cell  voltage  and  power  density  vs.  current  density  with  Ns  =  10.  (a) 
Dc  =  1  x  10-6  m;  (b)  Dc  =  lx  lO”5  m;  (c)  Dc  =  lx  10“4  m.  (*  )  T=  923,  (O)  T=  973  K, 
(>)  T=  1023  K. 


entropy  change  in  the  cell  reaction  under  normal  operating  con¬ 
dition  [22].  This  is  a  unique  thermodynamic  predominance  for 
DCFC. 

Different  from  gaseous  or  liquid  fuels,  the  concentration  of  car¬ 
bon  is  invariant  throughout  the  reaction  process  at  the  anode,  and 
the  phenomenon  of  fuel  crossover  is  inexistent  for  DCFC.  Therefore, 
Yf  is  equal  to  100%.  Flowever,  the  voltage  efficiency  yv,  is  affected 
by  the  operating  temperature,  anode  dimension,  particle  size,  and 
carbon  species,  etc.  Based  on  the  model  formulated  in  this  work, 
an  approximate  efficiency  of  DCFC  can  be  given.  When  the  current 
density  is  <200  Am-2,  the  overall  efficiency  y  is  higher  than  55%, 
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and  is  about  45-50%  when  the  current  density  is  >400 Am-2  in 
923-1023  K.  Finally,  it  should  be  noted  that  this  result  is  based  on 
the  DCFC  single  cell  using  graphite  as  fuel,  molten  Li/K  carbonates 
as  the  electrolyte  and  a  porous  Ni  cathode  as  the  one  used  in  a 
MCFC  [8]. 

4.  Conclusions 

The  performance  of  a  DCFC  was  simulated.  The  simulation 
results  reveal  several  important  feathers  of  the  very  promising 
DCFC. 

The  cell  is  under  the  control  of  anode  activation  polarization  at 
923-1023 1<  due  to  the  low  activity  of  graphite  used  in  this  work.  The 
working  temperature  has  a  decisive  effect  on  the  various  polariza¬ 
tion  losses,  especially  for  the  anode  activation  polarization,  which 
indicates  that  the  anode  surface  reaction  rate  is  a  key  factor.  The  per¬ 
formance  of  the  cell  depends  strongly  on  the  height  of  the  packed 
bed  anode,  and  the  carbon  particle  size.  This  means  that  the  reac¬ 
tion  surface  area  and  the  transportation  path  are  among  the  most 
important  factors. 

The  power  density  can  be  as  high  as  200-500  Wm-2  with  car¬ 
bon  particle  size  in  the  range  of  1.0  x  10-7  to  1.0  x  10-4  m  and  in 
923-1023  K.  The  approximate  total  efficiency  of  the  single  cell  is 
higher  than  55%  for  low  current  density  and  45-50%  for  high  current 
density. 
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